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Abstract. We present mid-infrared spectro-imaging (5 - 16/jm) observations of the infrared luminous interacting system Arp 
299 (= Mrk 171 = IC 694+NGC 3690) obtained with the ISOCAM instrument aboard ISO. Our observations show that nearly 
40% of the total emission at 7 and 15yum is diffuse, originating from the interacting disks of the galaxies. Moreover, they 
indicate the presence of large amounts of hot dust in the main infrared sources of the system and large extinctions toward the 
nuclei. While the observed spectra have an overall similar shape, mainly composed of Unidentified Infrared Bands (UIB) in the 
short wavelength domain, a strong continuum at ~ 13 fim and a deep silicate absorption band at 10 yum, their differences reveal 
the varying physical conditions of each component. For each source, the spectral energy distribution (SED) can be reproduced 
by a linear combination of a UIB "canonical" spectral template and a hot dust continuum due to a 230 - 300 K black body, after 
independently applying an extinction correction to both of them. We find that the UIB extinction does not vary much throughout 
the system (Ay < 5 mag) suggesting that most UIBs originate from less enshrouded regions. IC694 appears to dominate the 
infrared emission of the system and our observations support the interpretation of a deeply embedded nuclear starburst located 
behind an absorption of about 40 magnitudes. The central region of NGC 3690 displays a hard radiation field characterized by 
a [Nem]/[Neii] ratio > 1.8. It also hosts a strong continuum from 5 to 16/jm which can be explained as thermal emission from 
a deeply embedded (A v ~ 60 mag) compact source, consistent with the mid-infrared signature of an active galactic nucleus 
(AGN), and in agreement with recent X-ray findings. 

Key words. Stars: formation - Galaxies: individual: Arp 299 - Galaxies: individual: Mrk 171 - Galaxies: interactions - 
Galaxies: starburst - Infrared: ISM 



' 1. Introduction young, as shown by th e prominent 180kpc (130 long tidal tail 

h> ! (Iffibbard & Y unll999h . 

j_j ■ At a distance of 41Mpc (v he i = 3080 km s~', H = 

C3 ; 75kms _1 Mpc" 1 ), Arp 299 (Mrkl71 = IC694 + NGC3690) Arp 299 has been studied extensively at all w avelengths 

is one of the nearest interacting galaxies. Because of its prox- over t he past thirty years (see the original work bv lWeedmar] 

imity, its spatial extent (8 kpc), and its high infrared (IR) lumi- 1 19721). Most no t ably th ough, the mid-infrared and radio work 

nosity 1 (L IR = 5.16 x 10 11 L ), it is one of the prime candidates of Geh rz et aD Jl983|). the observations in the opti cal by 

for exploring the effects of triggered star formation activity as Augard e & Leaueuxl (I1985h. in t he ne ar-infrared by iGallaisl 

the tidal forces during an interaction lead to instabilities in the dl99lh and IWvnn- Williams et all Jl99ll) . al ong with the infor- 

galactic disks and rapidly funnel lar ge quantities of gas into the mati on on its mo l ecular gas kinematics by lAalto et all dl997l) 

dynamical centers of the galaxies JMihos & Hernauisllll996l and ICasoli et all i 19991) have provided ample evidence for 

ISanders & Mirabelll 19961) . Similar to the prototypical interact- the extreme characteristics of this system. More recently, the 

ing pair of galaxies NGC 4038/39, this system is dynamically high spatial resolution images ob tained with the Hubble Spac e 

Telescope (HST) in the near-IR jAlonso-Herrero etal1F2 000). 

as well as in the mid-IR using Keck JSoiferetalJl200ll) have 

Send offprint requests to: P. Gallais, pascal.gallais@cea.fr , , . . ,., , . . , . 

. , , . , T ' r f ... revealed numerous point-like sources that were undetected in 

Based on observations with ISO, an ES A pro ect with instruments , , ,• r , > 

f j j u cc a « u o. . / ■ ii .u tit . • c the optical, and firmed up our understanding of the complex 

funded by ESA Member States (especially the PI countries: France, r r i 



Germany, the Netherlands and the United Kingdom) with the partici- starforming history of the galaxy (see | Alonso-Herrero et al 

pation of IS AS and NASA 2000, an d references th erein). Following the notation suggested 

1 We use the standard definition of L IR (8 - 1000 ^m) = bvlGehrz et all \ 1 9831). the system is usually described bv the 

5.62 x 10 5 A 2 , r (13.48/ i2 + 5.16/ 25 + 2.58/ 60 + /ioo)Aa ( see g a l ax Y IC694 (source A) to the east, the nuclear region of 

ISanders & Mirabell 19961) . NGC 3690 (source B) to the south-west and two regions of 
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Fig. 1. a) The HST/NICMOS 2.2 //m image of Arp 299 from lAlonso-Herrero et aD fcOOOl) with the overlay of the ISOCAM 7 fim 
emission. The different components of the system are indicated. The 9 contour levels are set with logarithmic spacing between 1 
and 33 mJy arcsec~ 2 . The image scale is displayed with a bar of 2 kpc (10") at the top. b) Same as in a) but using the ISO CAM 
15 fim image as an overlay, having set the contour limits to 6 and 60 mJy arcsec~ 2 . Adapted from lCharmandaris et alJ -20021) . 



emission (sources C and C) to the north-west. Source B in 
NGC 3690 has been further resolved into several components, 
among which B 1 and B2 clearly dominate the global emission 
in the near- and mid-IR (see Fig. [T] for details). I nterest ingly 
it was also recently shown by ICharmandaris et al.l (J2002) that 
despite its inconspicuous appearance in the optical, IC 694 con- 
tributes nearly half of the infrared (8 - 1000 fim) luminosity of 
the system. 

One of the recurrent questions though in all studies of 
Arp 299 is whether t his ga laxy is a pure starburst as was pro- 
posed by WeedmanJ (119721) or whether some luminosity could 
be contributed by Active Galactic Nucleus (AGN). The fact that 
it is one of the most X-ray luminous galaxies (Ln.7_4.o kRV ~ 
10 42 erg s -1 , see lFabbiano et dll992_l_-ezas et al.lll998l) along 
with the flatness of its radio continuum emission had given 
some indications that an enshrouded AGN could be lurking 
in one of the nuclei. Only recently thoug h the hard X-ray 
(10 - 40keV) Beppo-SAX observations of iDella Ceca et alJ 
(2002) have provided unambiguous evidence that a deeply 
buried (JV H = 2.5 x 10 24 ctrT 2 ) AGN is located in one of the 
interacting components of the system. Furthermore, the same 
authors also demonstrated that the intrinsic AGN luminosity 
in both the UV and X-rays is almost an order of magnitude 
less than the infrar ed emission of the system calculated by 
IChar mandaris et al.l ll2002ll . indicating that the bulk of the in- 
frared emission of the galax y is indeed due to massive star for- 
mation as proposed bv Laurent et alJ (I20001) . The poor spatial 
resolution of the Beppo-SAX can not pinpoint the exact loca- 
tion of the AGN. Howev er, preliminary analysis of the Chandra 
data dZezas et alJl2003t) . which provide sub-arcsecond resolu- 



tion, suggests that the AGN might be located within source B 1 
of NGC 3690. 

Motivated by this recent activ ity, we decided to a nalyze 
in detail our 5-16 fim ISOCAM dCesarskv et al.lll996ri spec- 
trophotometric data on this system. Our observations provide 
good spatial resolution with unprecedent sensitivity and en- 
able us to examine in detail the dust properties of the obscured 
nuclear regions of Arp 299. The observations and data reduc- 
tion methods are described in Sect. 2, a brief description of the 
global mid-IR morphology of the system is presented in Sect. 3 
and the analysis of the mid-IR spectra is shown in Sect. 4. In 
Sect. 5, we discuss the implications of our findings on the phys- 
ical characteristics of the various infrared sources and our con- 
clusions are summarized in Sect. 6. 

2. Observations and Data Reduction 

Arp 299 was part of the ISO guaranteed time program 
CAMACTIV (PI. I.F. Mirabel) which had, as a prime goal, 
the study of the mid-IR propertie s of more th an 20 nearby ac- 
tive/interacting galaxies .Laurent et all koOO). It has been ob- 
served with ISOCAM with the Continuously Variable Filter 
(CVF), resulting in a full coverage of the spectral range from 
5 to 16 fim with a spectral resolution R ~ 40. The pixel size 
was 1.5" giving a total field of view of 48" x 48". The sys- 
tem was then fully mapped with an effective spatial resolution 
ranging from 3 to 4.5". The data reduction and analysis were 
performed using the CAM Interactive Software (CIA) 2 . Dark 

2 The ISOCAM data presented in this paper were analyzed using 
"CIA", a joint development by the ESA Astrophysics Division and the 
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subtraction was done using a model of the secular evolution 
of ISOCAM's dark current. Cosmic rays were removed using 
a multi-resolution median filtering method while the memory 
effects of the detector were corrected using the so-called IAS 
transient correction algorithm which is based on an inversion 
method ( Aber gel et alJl9 9&). The flat field correction was per- 
formed using the library of calibration data. Finally, individ- 
ual exposures were combined using shift techniques in order to 
correct the effect of jittering due to the satellite motions (am- 
plitude ~ 1"). These methods and th eir consequences are dis- 
cussed in detail in LStarck et al.l lfl999). We estimate that the un- 
certainty of our mid-IR photometric measurements is ~ 20 %. 

Experience with past ISOCAM observations has indicated 
that the overall pointing of ISO in the CVF mode is reli- 
able. However, given the disturbed morphology of the system, 
we verified the relative astrometry by comparing and iden- 
tifying the main features of our mid-IR images with those 
detected at other wavelengths after taking into account the 
shape of their spectral energy distribution (SED) as well as 
possible effects of absorption. For this purpose, we made 
an extensive use of the HS T WFPC2 and NICMOS data of 
lAlonso-Herrero et al,| (|2.QQ 0h- as wel l as the Keck 12 /mi im- 
ages obtained bv ISoifer et all J20Qll) . Although these mid-IR 
data are not as deep as our ISOCAM observations, their much 
higher spatial resolution clearly facilitated our efforts. Since 
one can not blindly rely on the astrometry information pro- 
vided by the ISO telemetry, we based our identification on the 
nucleus of IC 694 (source A), which is sufficiently isolated and 
its position is well defined in both the HST near-IR images as 
well as in the Keck 12 /mi maps. Finally the stability of the 
roll angle of the telescope during the observations and the apri- 
ori knowledge of the separation of the brighter mid-IR point 
sources enabled us to securely identify the position of the re- 
maining components. 

3. The Mid-IR morphology 

As mentioned earlier, the compar ison between our data and the 
high resolution Keck images of lSoifer et alJ J200ll) allowed us 
to clearly locate the four components A, B, C and C as the 
dominant sources of emission in our ISOCAM images. These 
identifications further revealed that the centroid of our source B 
actually corresponds the position of B 1 . Even though the angu- 
lar separation between Bl and B2 is only ~ 2", well below our 
spatial resoluti on, the fact that B2 is only marginally detected 
in the maps of ISoifer et al leads us to believe that its 

contribution to the mid-IR emission of B is minimal longward 
of 7 /mi. Thus the conclusions we will draw for source B actu- 
ally reflect the characteristics of B 1 . 

In an attempt to increase our sensitivity to faint emission 
features, we used the CVF dataset to construct broad band im- 
ages at the two most commonly used ISOCAM filters centered 
at 7 and 1 5 /mi. The resulting 5-8.5 /mi and 12-18 /mi images, 
which correspond to the LW2 and LW3 filters, reach a lcr-rms 
noise of 0.27 and 0.44 mJy arcsec~ 2 respectively. They are pre- 

ISOCAM Consortium led by the ISOCAM PI, C. Cesarsky, Direction 
des Sciences de la Matiere, C.E.A., France. 



sented in Fig.^ adapted from lCharmandaris et alJ (J2002), with 
their contours overlaid on th e NICMOS HST data obtained by 
lAlonso-Herrero et alJ |2000 ) . The LW2 image is dominated by 
the contribution of the Unidentified Infrared Band (UIB) fea- 
tures and as a result traces mostly the contribution of photo- 
dissociation and quiescent star-forming regions. The LW3 im- 
age, on the other hand, mainly samples the thermal signature 
of hot dust from massive star-forming regions as well as emis- 
sion from high excitation forbidden lines such as [Nen] 12.8 /mi 
and [Nem] 15.6 /mi. A more detailed description on how the un- 
derlying mid-IR SED affect s the LW2 and LW3 colors in ac- 
tive galaxies is presented in lLaurent et al ] J2000l) . The mid-IR 
fluxes of the various sources are reported in Tabled 

In addition to the four main components (A, B, C and C) 
detected in our mid-IR images, we note that there is consid- 
erable diffuse emission, in particular at 7 /mi, in the region be- 
tween IC 694 and NGC 3690 as well as around the nuclei of the 
galaxies. Our equivalent broad band images reveal that this dif- 
fuse component, which was undetected in the shallower ground 
based images, accounts for 44 % and 36 % of the total 7 and 
15/mi fluxes originating from Arp299. This diffuse emission 
is shown in Fig.|2]which displays normalized profiles at 7. 1 /mi 
(continuum) and 7.7 /mi (UIB band) of sources A and Bl com- 
pared to the PSF profiles. Together with the profiles, images 
of the diffuse emission at these wavelengths are presented, ob- 
tained by subtraction of a synthetic image made of PSFs posi- 
tioned at the locations of sources A, Bl, C and C to the images 
of the system in the 7.7//m UIB and in its continuum. If we 
derive an equivalent IRAS 12 /mi flux assuming a power-law 
SED in the 5-18 /mi range, we find that 60 % of the overall 
IRAS 12 /mi flux is accounted for by the four sources, with A 
and B 1 being the dominant ones (22 and 28 % of the total flux 
respectively). 



Table 1. Broad-band mid-IR photometry of Arp299 



Source 


LW2 1 (mJy) 


LW3 1 ' (mJy) 


LW3/LW2 


A 


325 


1860 


5.72 


B 


505 


1951 


3.86 


C 


126 


461 


3.66 


C 


76 


232 


3.05 


Total 


1846 


7037 


3.81 



f using a beam of 4.5" in diameter (aperture corrections were applied 
to account for the overall extension of the PSF). 

All main components in the system are very bright in both 
filters and have extremely red colors. The 15 to 7 /mi flux ra- 
tios LW3/LW2 range from 3 to 5.7 and they are quite high 
compared with what is typically found in normal g alaxies 
where they vary from 0.7 to 1.2 jRoussel et alJl200ll) . Even 
in the case of the Antennae galaxies, Knot A, the brig htest 
15 jUm region, exhibits a LW3/LW2 ratio of 2.6 ( Mi rabel et alJ 
Il998t) . while the same indicator in the massive star-forming 
complex within the outer r ing of the Cartwheel galaxy is 5.2 
JChar mandari s et alJ Il999h . In the absence of an AGN-type 
activity, such high ratios in Arp299 indi cate a star-forming 
mechanism of a remarkably high efficiency ( Gall ais et all 1 999: 





Fig. 2. Normalized surface brigthness profiles of sources A (panels la, lb) and Bl (panels 2a, 2b) in the continuum at 7.1 fim 
(dot-dashed line with diamonds) and in the UIB band at 7.7yum (filled dots on continuous line) overplotted to the normalized 
profile of the PSF (dashed line) along the 2 detector axis (considered pixels are tagged by a line in image 3b). Image 3a and 
3b represent synthetic image built with adjusted PSFs positioned at the locations of A, Bl, C and C (3a: full dynamical range; 
3b: 5% of the maximum). Images 3c and 3d display the extended emission respectively in the 7.7 /^m UIB and in the 7.1 /urn 
continuum, i.e. the residual emission after subtraction of the synthetic image of PSFs to the image of the system in both bands 
(same dynamical range than the one applied for image 3b) 



Laurent et all200(31) . a result which has a lso been suggested by 
Lai et alJ i 19991) and iRouan etall i 19991) based on their near- 
IR adaptive optics observations. Using a gray body model to 
fit the mid- and far-IR SED of the svstem. lCharmandaris et alJ 
(2002) have estimated that the infrared luminosities of the three 
main components (A, B1+B2 and C+C) range from 0.44 to 
1.8 x 10 11 L Q . This implies that each of these sources is several 
times more luminous than most starburst galaxies observed in 
the local Universe such as the prototypical M 82. 

Examining the individual sources in more detail, it is strik- 
ing that, within the limits of our resolution, component A is 
dominated in the mid-IR by a strong nuclear source with a faint 
extension towards the southeast. This nuclear source is com- 
pletely absent from the optical HST ima ges and only becomes 
appar ent in the original near-IR maps of Wvnn-Will iams et alJ 
lll99ll) . In the near-IR, this is interpreted as the emission from 
the old stellar population in the bulge of IC 694. However, the 
mid-IR spectrum (see Fig.|4ji differs clearly from the one of a 
population of old stars, such as in ellipti cal galaxy, in w hich 
there is no emission or absorption band llAfhev et alJl2002l) . It 
is rather similar to th e one we encounter in star forming regions 
( Laurent et alJ fcoOO). In the mid-IR, the dominant process is 
star formation which is not evidenced at shorter wavelengths 
(visible and near-IR) and it appears mainly concentrated in 
the nuclear region. At longer wavelengths, the emission of A 
dominates the global infrared emission of the system, mak- 



ing A the most lum inous component of the system at 38 /im 
JCharmandaris et aTll2002t) . This implies that even though star 
formation does happen over an extended area of IC 694 - as ev- 
idenced by the dust lanes and star forming tidal s treame rs vis- 
ible in the HST optical imagery of Malkan ~et alJ Jl998l) - and 
does contribute in the mid-IR by creating an extended emission 
seen in Fig.^ the bulk of the infrared flux originates from its 
nucleus. 

The mid-IR structu re of NGC 3690 is more com plex. As 
discussed in detail by lAlonso-Herrero et alJ (|2000), a large 
number of point-like sources can be seen in the near-IR, 
with Bl and B2 largely dominating the overall emission of 
the galaxy at these wavelengths. The spectral properties of 
these two main components though present a remarkable con- 
trast. As the wavelength increases, Bl becomes progressively 
brighter while the contribution of B2 - clearly the most lumi- 
nous s ource of NGC 3690 in the visible and in the J and H 
bands ( lAlonso-Herrero et all200ol) - gently declines longward 
of 1.6jum. Furthermore, the latter is bare ly detecta ble in the 
high resolution mid-IR images o flSoifer et alJ 11200 II) . while Bl 
clearly unveils a very luminous core at 12//m. The red colors 
of Bl and its d ifferential extinction of Ay = 15 mag relative 
to B2 lAlonso-Herrero et alJl200oT) are indicative of its much 
more obscured nature and reveal the presence of large concen- 
trations of dust and gas in this area. As a side note, it is worth 
mentioning that Bl is completely invisible in the HST/FOC 
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220 nm UV map of the galaxy (see F ig.E), and it is very faint 
at 300 nm (see Windhorst et al. 2002Fig. 3.34a,b), an additional 
evidence of its high dust content. 



N 2" = 400 pc 




Bl 

NGC 3690 



Fig. 3. UV image (220 nm) o f NGC 3690 taken w ith the FOC 
on-board the HST (see also iMeurer et alJ ill 9951) 1 with posi- 
tions of B 1 and B2 marked. One can see that no emission at 
this wavelength is associated with B 1 confirming its highly ob- 
scured nature. 

Sources C and C are clearly resolved at 7 //m but the emis- 
sion from C diminishes relatively to C at longer wavelengths 
(see Fig.0. It has been proposed that source C could be the 
nucleus of a third galaxy in the interacting system, but this hy- 
pothesis is questioned by the l ocation of this source on the [J- 
H]/[H-K] color-color diagram (Gallais 199JJ), far from the lo- 
cus of the other nuclei. Recent CO observations jAalto et alJ 
ll997HCasoli et alJ l999) also suggest that C and C are dynam- 
ically linked and can be considered as star-forming complexes 
in the overlap region of NGC 3690 and IC 694. 

4. Mid-IR spectral properties and absorption 

In Fig.|4] we present the 5 - 16 /an spectra of the four main 
sources A, Bl, C and C of Arp299. They have been obtained 
with a 4.5" x 4.5" aperture centered on the brightest spot of 
each component. A wavelength-dependent scaling factor has 
also been applied to correct for the extended flux of the PSF 
lying outside of the aperture used. We observe that the overall 
shape of each spectrum is similar, with the UIB features dom- 
inating the emission in the short wavelength range and with 



a warm continuum becoming visible longward of 10 /an. Our 
measurements are also consistent with the 8 — 13 /an spectro- 
scopic observations already obtained by Dudlev ( 1999). 

A more thorough analysis though, reveals distinct underly- 
ing physical properties of each region. While the mid-IR SED 
of source C is typi cal of what is gen erally observed in star- 
burst galaxies (e.g., Laur ent et al . 2000), the spectra of sources 
A and C display very strong absorption in the 9.7 /mi silicate 
band and a significant contribution of the thermal continuum at 
A > 12 /an. Sources C and C also display a strong [Nem] emis- 
sion line, which as we discuss in section 5, is not as apparent 
in the other components. Finally, we note that the spectrum of 
source Bl is rather intriguing since it is dominated by an im- 
portant hot dust continuum throughout the whole CVF spectral 
range, with only a minor contribution of the UIB features. 

To provide a more quantitative insight into the properties 
of these sources, their 5 - 16//m SEDs were carefully ana- 
lyzed and compared with a library of simulated-CVF spectra. 
The latter were created with a two-component model described 
hereafter. The role of the first component is to take into ac- 
count the presence of the UIB features, usually exhibiting a 
canonical SE D as it has been observed in vari ous star-forming 
galaxies JPale et all200lURoussel et all200ll) . so we used the 
average 5-16 //m SED of the M 5 1 galactic disk as a template 
spectrum. The second component, which reproduces the warm 
continuum observed in the long wavelength part of the CVF, 
was modeled using a black body with a temperature typically 
ranging from 150 to 300 K. We also considered the use of a 
power law with a spectral index varying from a — 1 to a — 2, 
as well as black bodies modified by an eirrissivit y law (see 
also lSturm et alJ2000HForster Schreiber et al.l20 03). However, 
given the similar shape of the power law and black body func- 
tions at these wavelengths and temperature regimes as well as 
the limited wavelength coverage of our spectra, we found that 
we were unable to discriminate between one of the two partic- 
ular spectral shapes. As a result, we will restrict our analysis 
using simple black body fits. Selective reddening was then ap- 
plied independently towards thes e two comp onents (similarly 
to the approach followed by Tran et al. 200 11), and a variety of 
extinction laws dMathisI 19901: [Dudley & Wvnn- Williams! 19971 
iLutzll 19991 iLi & Drainefe OO 1 ) were explored assuming either 
a uniformed dust screen geometry (F Aou , = Fx,m X exp~ Tl , 
ta = A^/1.086) or a homogeneous mixture of dust grains and 
sources (F /lml , = F AJ „ x (1 - exp -1 "- 1 )/^ , t a = A^/1.086). 
Varying the fraction of the UIB and continuum emission to 
the total mid-IR flux, a collection of CVF spectra was subse- 
quently created and ^-tested against our data to select the best 
simulated SED available in our library. Note that the strong 
ionic emission fine s usually observed in starburst environments 
dSturmet all EoOO) such as [Nen] at 12.8 /an and [Nem] at 
15.6/an were not taken into account in our simulated CVF 
spectra as they are not included in our model. Moreover, at 
the spectral resolution of our CVF, [Nen] is blended with the 
12.7 /an UIB feature. Therefore the 12-13.2 //m and 15-16 /mi 
wavelength ranges were omitted in our^ 2 calculations. 

In Fig.|4] we present the best fits we obtained for the four 
ISOCAM spectra, along with the decomposition into their two 
components. While the specific results which were derived for 
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Wavelength [wm] Wavelength [|xm] 




Wavelength [/im] Wavelength [/4m] 

Fig. 4. Rest-frame mid-infrared spectra of the four main sources (A, Bl, C and C) in Arp299 (solid lines). They were obtained 
with a 4.5" x 4.5" aperture centered on each component, and they were corrected via a scaling factor to account for the extended 
flux of the PSF lying outside of the aperture. For each component, we have also plotted our best fit (dash ed line ) obtained w ith the 
combi nation of a UIB template and a black body continuum (dotted lines ) extinguished with the law of Dudl ey & Wvnn- Wi lliams 
On each panel, we indicate the physical properties of the continuum: its contribution to the mid-IR flux, the absorption 
magnitudes followed by the dust obscuration geometry (screen against mixed model), and the temperature of the fitting black 
body. We also found that the extinction level of the UIB component is negligible (Ay < 5 mag) in the four sources. Note that the 
[Nen] 12.8 /im and [Nein] 15.6 /mi emission lines were not taken into account in our models. 



each source will be more thoroughly discussed in the next sec- 
tion, we highlight hereafter some general points of our findings, 
as well as a number of limitations related to our fitting proce- 
dures. 

One of our most striking results is that our best fits indicate 
an elevated extinction for the underlying continuum emission 
of sources A, B 1 and C and a very low extinction (A y < 5 mag) 
of the template which contributes to the UIB feature emis- 
sion. This may suggest not only a physically different origin 
but also a different spatial distribution for the two components 
observed in our spectra. Such a picture could actually be un- 
derstood if we consider that the continuum emission origi- 
nates from deeply dust-enshrouded regions while the UIBs are 
more diffuse and have a widespread and surface distribution. 
Since the UIBs mostly trace areas of photo-dissociation sur- 
rounding Hn regions, this would not be so surprising given 



the filamentary optical/near -IR emission (Mal kan etafl ll998: 
Alonso-Herre ro et alj |2000 |) and the extend ed gas streamers 
(Sargent & Scovill eTl99lUCasoli et alJll999l) which have been 
observed in the in-between regions of Arp 299. 

A potential limitation in our approach was revealed by the 
degeneracy originating from the unknown geometry of the dust 
in the obscured regions. For all sources but Bl, we found that 
both the screen and the mixed models can actually lead to sim- 
ilarly good overall fits to the observed spectrum in terms of 
the x 1 minimization, even though the underlying assumptions 
for the geometry of the sources are different. For example, the 
continuum emission of source A can easily be reproduced ei- 
ther with a 230 K black body extinguished with Ay = 40 mag 
in a screen geometry (see Fig.|4}, or with a much cooler com- 
ponent at T = 160 K with Ay = 140 mag assuming the mixed 
model. However, we wish to stress that the physical meaning 
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of such single-temperature fits is rather uncertain since i) the 
observed spectrum is the result of no single grain population 
but a linear superposition of grains with a specific size distri- 
bution (see lDale et all200ll) . ii) the dust grains responsible for 
this warm continuum are not in thermal equilibrium, and Hi) 
the presence of temperature gradients in the close vicinity of 
the mid-IR sources is quite likely. Therefore, a more meaning- 
ful approach would rather be to compare the four fitted spectra 
to each other after having primarily defined a given geome- 
try for the dust obscuration (see Sect.[5]for such comparisons). 
Furthermore, we note that both the fraction of UIBs to the to- 
tal mid-IR flux and their negligible extinction measured using 
our fitting procedure are roughly independent of the screen or 
mixed model assumption. 

Table 2. Intensities and ratio of the [Nen]12.8/mi and 
rN eml 15.6 um for Arp299. The values of [Nen] were measured 
bv lDudlevi (119991) and should be considered as upper limits as 
they may be contaminated by the 12.7 //m PAH feature (see 
text). 



Source 


/([Nen]) 

(KT l5 WnT 2 ) 


/([Nem]) 

(1(T 15 WitT 2 ) 


/([Nelll]) 
/([Nell]) 


A 


< 2 






Bl 


0.7 


> 1.3 


> 1.8 


C 


2.0 


1.3 


>0.65 


c 


< 0.9 


0.5 


>0.55 



To investigate the hardness of the radiation field in the dif- 
ferent regions of the system we can use the ratio of the ionic 
Neon lines [Nen]12.8/im and [Nem]15.6yum. Even though we 
can easily measure the [Nem] line, the spectral resolution of the 
CVF (R ~ 40) is just too low to distinguish the [Nen] 12.8 yum 
line from the 12.7//m PAH band. As a result we decided 
to use the [Nen] values of the system measured by iDudlevI 
(1999) who calculated the intensity of the [Nen] lines in all 
sources of Arp299. His R ~ 60 spectrograph can almost re- 
solve the [Nen] line, and the 5.5" circular aperture he used 
is only ~ 25 % smaller than the aperture used to extract our 
ISOCAM/CVF spectra presented in Fig.|4] Our results are pre- 
sented in Table |3 In principle, the values of [Nen] may still 
contain some contribution from the 12.7 fim PAH feature, thus 
making the [Nem]/[Nen] ratio a lower limit. However, this con- 
tribution must be small since the sum of line intensities mea- 
sured by Dudlev ( 1999) in the regions Bl, C, and C is con- 
sistent with the high spectral resolution ISO/SWS measure- 
ments of Tho rnlev et alJ J2000h . We note an inconsistency be- 
tween the sum of our measurements o n B 1 , C and C and th e 
integrated measurement presented in iThornlev et ail J2000t) . 
Whether their value underestimated or our overestimated is not 
easy to determine, and spectroscopic observations of this target 
with SIRTF will resolve this issue. 

One of the great advantages of the ISOCAM/CVF observa- 
tions is that we can extract narrower-band images (AA ~ 1 fim) 
of our target, effectively constructing two dimensional maps 
of specific spectral features. Several maps of selected spectral 
regions were created for Arp 299 and are displayed in Fig.|5] 



The first is the hot dust continuum at 5 fim, which could be 
due to either thermal emission from an AGN or the Rayleigh- 
Jeans tail of the photospheric emission from an old stellar pop- 
ulatio n, or fluctuating s mall dust grains without aromatic fea- 
tures ( He lou et alJ2 000). The integrated emission of the 6.2 fim 
UIB is also presented, since it isolates the emission from the 
photo-dissociation regions of the merger and despite the fact 
that it is weaker than the 7.7 fim UIB it is much less affected 
by the depth of th e 9.7yum silicate band (see discussion in 
lLaurent et alJ l2000). We also mapped the total flux within the 
band of the silicate absorption at 9.7/mi, and that of the ionic 
[Nem] 15.6 /mi emission line which is indicative of the hardness 
of the radiation field (e.g.. lThornlev et afll2000l) . For the latter 
and the 6.2 fim UIB, we identified the corresponding spectral 
feature within each pixel of the CVF images, and carefully re- 
moved the underlying continuum so as to only keep the con- 
tribu tion of the given em ission line ("see lLaurent et alJ [2000 
and lLeFloc'hetalJbOOll for an illustration of the method). 
The dominant sources of mid-IR emission (A, B and C) are 
clearly apparent in all maps, even though their relative strength 
varies at each wavelength. However, the image of the ionic 
[Nem] 15. 6 /mi line remains indicative of the spatial distribu- 
tion of this feature. The spectra of individual pixels are noisier 
than any integrated spectra and the determination of the con- 
tinuum pixel-per-pixel can be altered by this noise. Moreover, 
especially in the case of source Bl, the depression in the SED 
over 16 fim can lead to a bad estimation of the continuum and 
then to a bad determination of the flux in the line. So, inte- 
grated spectra on each main source are more reliable in term 
of relative flux. Even if C appears as the dominant source in 
the [Nem] 15. 6 //m image of Fig.[5]d, Bl remains the brighest 
source in this line. 



5. Discussion 

5.1. The nuclear region of IC 694 

In our mid-IR images, IC694 (source A) appears as a 
bright source corresponding to the r adio peak observed at 
6 and 20cm by |Gehrz et all dl983l). As pointed o ut (see 
lAlonso-Herrero et alJl2000t ICharmandaris et aDl2002l and ref- 
erences therein), unlike the infrared view, images in the UV 
and visible do not reveal any clear optical counterpart at this 
position, while large quantities of high de nsity molecular gas 
are concentrated within the central 500 pc (Sar gent & Scovillel 
ll99lUCasoli etaljri 999). This gas fuels a massive starburst as 
evidenced by a supernova rate of 0.65 yr~', nearly five times 
the c ombined rate observed in B and C ( lAlonso-Herrero et alJ 
2000). This large gas reservoir can easily explain the high ex- 
tinction toward the nucleus. In the optical, a dust lane runs in 
the southeast/northwest direction and "patchy" emission from 
surface star formin g regions is se en in an extended area around 
the nucleus (Malk an et alJll998l) . Not surprisingly, those same 
regions contribute to the 6.2 fim UIB emission and, as we ob- 
serve in Fig.[5p, the angular size of source A is considerably 
larger at this wavelength (see also Fig.|2]which shows this dif- 
fuse emission). 
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Fig. 5. Inverted grey-scale images of Arp 299 in selected spectral regions, depicting a) the hot continuum emission observed at 
5 /mi, b) the photo-dissociation regions obtained by integrating the mid-IR emission within the 6.2/mi UIB, c) the continuum 
emission integrated over a AA ~ 1 //m bandwidth in the range of the silicate absorption at 9.7/rni, and d) the distribution of 
the ionic [Nem] emission line at 15. 6 /mi, after subtraction of the underlying continuum. This last image must be considered 
carefully because the continuum per pixel on source B may have been badly estimated due to the low signal to noise ratio of the 
emission in the line and the probable influence of the wide 18 /mi silicate absorption band. For each image, the overlaid contours 
have been set with logarithmic spacing. By increasing order of wavelengths, the minimum levels are respectively 0.44, 0.18, 0.22 
and 0.49mJy arcsec~ 2 while the maxima are 10.7, 2.0, 15.5 and 1.5 mJy arcsec~ 2 . The positions of the brightest sources (A, Bl 
and C) over the 5-16 //m wavelength range are marked with crosses. The four images have the same orientation. The average 
FWHM of the ISOCAM PSF over the full CVF spectral range (~ 4") is indicated in panel a) as well as the scale in arcsec and 
kpc. 



Our mid-IR spectrum of this source (see Fig.@Ji reveals a other component of the system. Moreover, the 9.7 //m silicate 
strong continuum with a slope past 10 /mi steeper than in any band appears saturated. The continuum emission from A can 
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be fitted with a 230 K black body absorbed by Ay ~ 40 mag as- 
suming a uniform screen of dust. We also note that the longer 
wavelength part of the spectrum seems to decrease again past 
15.5 /mi. The same behavior can be seen in source B. Since 
we do know from multiple observations between 20 and 38 fim 
that the spectrum of A rises quickly maki ng it the dominant 
sourc e of the IR luminosity of the system JCharmandaris et al] 
120021) . it is reasonable to attribute this brief decrease in flux 
to the presence of the 1 8 um silicate feature. Indeed, as shown 
bv budlev & Wvnn-Williams1(ll997h . the two silicate bands are 
proportional in strength/depth so, given the observed shape of 
the 9.7 fim feature, the 18 /mi one would have to be strong too. 
Unfortunately given the limited wavelength coverage of the 
CVF and the fact that memory effects of the detector are more 
pronounced at the longer wavelengths which are observed at 
the beginning of the scan, we can not elaborate more on its ac- 
tual strength. We will have to wait for the Infrared Spectrograph 
(IRS) aboard the Space Infrared Telescope Facility (SIRTF) 
which will provide complete wavelength coverage up to 40 fim 
to address this issue more accurately. 

A unique aspect of IC 694 that separates it from all other 
areas in the galaxy is the weak [Nem]15.6//m emission. As 
we can see from Fig. |4^, even though a suspicious feature at 
15. 6 /mi could be the [Nem] line, its equivalent width is very 
small and makes it practically impossible to measure with our 
resolution. This can not be easily explained by advocating that 
the blue wing of the 18 /mi silicate feature extends to the loca- 
tion of [Nem] reducing the line flux. If simple extinction was 
the culprit then how could one detect the line in both regions B 1 
and C where even higher values of Ay are found? Furthermore, 
why does it seem that the spatial distribution of the faint [Nem] 
avoids the nucleus but is associated with the dusty filaments to 
the southeast of source A? 

The fact that the mid-IR neon lines of IC 694 have also 
been obs erved with IS O/SWS with superior spectral resolu- 
tion bv lThornlev et alJ J2000I) can help us explore these ques- 
tions. Using a considerably larger aperture (14 " x 27 ") these 
authors easily detected both [Nen]12.8//m and [Nem] 15. 6 fim 
toward IC694. They measured a [Nem]/[Nen] ratio of 0.29, 
a factor ~ 2.5 times lower than what was found for the 
whol e of NGC 369 (combined sources B and C). Modeling 

I 1 I rl 

by Thornlev et al. (2000) has shown that low values of the 
[Nem]/[Nen] ratio can be explained in terms of a larger age 
of the starburst, lower limit in the upper mass function cut off, 
higher metallicity, or a reduced ionization parameter U. Our 
observations are actually consistent with all these possibili- 
ties if we consider that IC 694 contains the larger fraction of 
molecular gas and cold dust and dominates the IR luminosity 
of Arp 299. This implies that the starburst in source A has been 
taking place for a longer period of time than in sources B and 
C since the former seems to be closer to dynamical equilib- 
rium than the two sources in NGC 3690 which ar e even more 
distur bed and obviously still strongly interacting dCasoli et alJ 
Il999i) . One could further speculate that the geometry in A i s 
different from the single cluster case of Thornle v et alJ (12000). 
The spatial extent associated with the age of the starburst indi- 
cates that the distribution of the enshrouded stars is more ran- 
dom forming clusters resembling the mix gas/star cluster case. 



Furthermore, more generations of stars must have been created 
in IC 694, and passing through their AGB phase, they produced 
the dust we observe and also lead to a higher metallicity in the 
region. As we will discuss in the following section this is not 
the case for C which mu st harbor the youngest OB stars in the 
system dSoiferet al.l200ll) . 

If the above scenario is correct then the reason why we de- 
tect [Nem] in the southeast outskirts of source A, and not in the 
center, is simply due to the fact that these are the regions wher e 
young stars form along the tidal filaments (Mal kan et al.ll9 98). 
In these filaments, the geometry is much more favorable so that 
"bubbles" in the interstellar medium form easier exposing the 
high ionization radiation of some underlying recently formed 
massive stars. 

5.2. On the mid-IR properties of sources C and C 

The shape of the mid-IR SED of these regions (see Fig.@J is 
consistent with those of star-forming regions in other extra- 
galactic sources, where a hot thermal continuum strongly af- 
fects the PAH emission. Furthermore, the values of the mid- 
IR diagnostic flux ratios LW3/LW2 (see Tabled are similar 
to th ose found in regions dominated by young massive stars 
(see [Laurent et alJl2000I) . As a result, our data provide fur- 
ther s upport to the conclusions of other authors dSoifer et alJ 
l200ll and references therein) that C and C host extraordinary 
extra-nuclear starbursts. 

Although it is probable that some emission from C could 
contribute to C because of their small angular separation, the 
fact that their measured absorptions are different indicates this 
is not a very likely event. Our modeling of the extinction sug- 
gests that the environments for both sources C and C are sim- 
ilar, with a higher extinction toward C. For both, as for source 
A, the observed spectrum can be reproduced by the simple su- 
perposition of a UIB contribution and a black body heated to 
~ 230 K in order to account for the hot dust continuum which 
dominates above 10 fim. The estimated extinction for sources 
C and C is Ay of 10 and 50 mag respectively assuming a fore- 
ground screen geometry. 

The hardness of the UV radiation field in these sources can 
be characterized by the [Nem]/[Neii] ratio. Source C presents 
a ratio of > 0.65 and C a ratio > 0.55 3 , bo th greater than 
the on e measured in component A (0.29 from iThornlev et alJ 
2000). This would indicate that the UV radiation field, and con- 
sequently the fraction of young massive stars is higher in C and 
C than in A, leading to younger starb ursts in these re gions, in 
agreement with the conclusions of Soif er et al . (2001 ) as well 
as of Alonso-Herrero et al. (2000) who derived ~ 5 Myr, 4Myr 
and ~ 1 1 Myr for C, C and A respectively (based on an evo- 
lutionary starburst model using parameters derived from their 
observations). 



3 We must be careful considering [Nem] in source C. We notice that 
the emission in this line is only marginally mapped at the C position 
(see Fig.|5J. Because C and C are separated only by 5", the feature 
we observe in the C spectrum may be a relic of the extended wings of 
the PSF (FWHM ~ 4.5") centered on source C. 
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5.3. The nature of source B: Is there mid-IR evidence 
for an AGN? 

As we mentioned earlier, Bl appears to be the brightest source 
of the system in the 5 to 16 /mi spectral range and remains unre- 
solved by our observations. While our spatial resolution is not 
sufficient to unambiguously resolve B 1 and B2, the fact that 
the centroid of the source located at B posi tion coincides very 
well with the 6 cm radio source observed bv lGehrz et alJ Jl983') 
indicates reasonably that we actually map mainly the emission 
coming from Bl. This is in agreement with higher spatial reso- 
lution ground-based observations from lSoifer et alJ 1 1200 ll) who 
resolved Bl and B2 but showed that, d espite the fact that B2 
dominates the flux below 2//m (see also lWvnn-Williams et alJ 
1 199 it lGallaislll99ll) . its emission decreases with increasing 
wavelength and it is actually ~ 35 times fainter than Bl at 
12.5 /mi. Consequently, all our conclusions on the mid-IR char- 
acteristics of source B=B1+B2 would actually reflect the prop- 
erties of B 1 . 

Until recently, most researchers had focussed on the star- 
burst properties of this region, given the absence of any evi- 
denc e for an active nucleu s. Based on the Ha ( see Fig. 3 of 
Al^jns^HeEgr^etalJ^OO) and CO observations JCasoli et alJ 
1999t lAalto et al J 1 1 997n 7it is clear that Bl contains a large 
amount of ionized gas, as well as cold molecular gas which 
could sustain massive star formation activity for more than 10 8 
years. Our observations though indicate that the physical char- 
acteristics of Bl, as displayed by its spectrum, make it stand 
out in the mid-IR in two ways. 

First, Bl displays the strongest continuum emission over 
the whole 5-16 //m range with only a weak contribution due 
to PAH emission (see Fig.|4}. Using our model, we find that 
this continuum can be best fitted with a 300 K black body ex- 
tinguished by Ay ~ 60 mag in a mixed-model geometry and 
accounts for 85 % of the total mid-IR emission of B 1 . Despite 
the high extinction derived from our model, Bl remains the 
brightest source of the system in the 9.7 //m silicate absorption 
band and its "hot continuum" at ~ 5 //m is 1 18 mJy, almost 5 
times higher than source A which dominates the global IR lu- 
minosity of the whole system. 

Second, the [Nem]/[Neii] ratio of this source is at least 
1.8 (see Table Et , a value cons iderably greater than the ratios 
measured in the Thornl ev et alJ (12000) starburst sample which 
typically ranges from 0.05 to 1, with the only exception of 
the low metallicity systems NGC5253 and IIZw40 (with val- 
ues of 3 and 12 respectively). Clearly, the radiative environ- 
ment in Bl deviates from the one seen in classical starburst 
galaxies, and it appears to resemble more wha t is observed in 
young galactic Hn regions such as W51 (Thorn lev et al.l2 000). 
ISturm et alJfe002l) reports measurements in the Neon lines for 
Seyfert galaxies which clearly indicate a higher ratio in pure 
Seyfert galaxies (usually > 1 with a maximum of 2.8 in their 
sample) than in mixed Seyfert/starburst systems (< 0.6 typi- 
cally). On the other hand, the ratio [Nem]/[Oiv] is usually low 
(< 1.5) in Seyfert gal axies as computed from values reported 
bv lSturmetatl ( l2002l) while it is hig her in starburst galaxies, 
ranging from 5 to > 30 Csee lVerma etall 12003). Assuming 
that C and C contributes to the flux measured in [Oiv] by 



IVerma etaD ( l2003l) . Bl then presents a ratio [Nera]/[Oiv] > 13. 
Even if those ratios are not criteria used to discriminate Seyfert 
and starburst galaxies, the association of a high [Nem]/[Neii] 
with a high [Nem]/[Oiv] in source Bl makes it of very intrigu- 
ing origin. 

The fact that no UV emission is detected from B 1 is fur- 
ther evidence of the high obscuration in the region. In this con- 
text, the PAH features could be attributed to individual heat- 
ing sources distributed around Bl, in a diffuse component on 
the line of sight, in such a way that the filling factor of their 
corresponding Stromgren spheres is high and hence there is 
a smaller volume available for th e photo-d issociation r egions 
to contribute to the PAH flux ("see lAlonso-Herrero et alJ feOOO. 
Fig. 7b). 

However, hard UV photons can also be produced by an 
AGN. Based on a diagnostic diagram which compares the ratio 
of the continuum at 14- 15 //m to the continuum at 5.1 -6.8 /mi 
versus the ratio of the streng th of the 6.2 jum feature to the con- 
tinuum at 5.1 -6.8 jum. lLaurent et alJ J2000b attempted to quan- 
tify the AGN/starburst contribution in the mid-IR. For source 
Bl, their Fig. 6 suggests that almost 40 % of the mid-IR emis- 
sion could originate from a very embedded AGN. The tell-tale 
sign would be the emission observed via the "hot continuum" 
at 5 /mi as illustrated in our Fig. [5] Even though the mid-IR 
evidence was not strong enough to propose the presence of an 
AGN detection, their suggestion was in fact proven true beyond 
any doubt with the recent Beppo-S ax X-ray observations of 
Arp299 by iDella Ceca et all J2002I) . These authors explained 
the hard X-ray emission and spectral shape/luminosity of the 
system with a buried Compton-thick AGN of intrinsic lumi- 
nosity Lo.5-iookeV - 1.9 x 10 43 erg s , completely absorbed in 
the visible and in the IR (Ay > 1 000 magnitudes). Their result, 
as well as fu rther w ork based on Chandra data which places the 
AGN in Bl JZezas et alJ2003l) . readily explains our inability to 
probe the center of B 1 in the optical and the UV. 

An interesting point which remains to be explored in the 
future is whether the AGN actually defines the location of a 
galactic nucleus in NGC 3690. The reason for this question is 
that B2 is the most extended source in the visible as seen in 
the HST WFPC2 (0.59 //m) images and the brightest one in the 
near-infrared, while Bl remains unresolved or marginally re- 
solved ( Alonso-He rrero et al . 2000). Moreover, unlike Bl, B2 
does display a high CO index characteristic of a population of 
late-type stars (Satv apal et aDf l999>. Based on its near-IR col- 
ors and its size at these wavelengths, B2 appears to be mainly 
composed of a mix of old stars, gas and dust. If this source is a 
relic of an older single isolated starburst region outside the nu- 
cleus traced by B 1 , it should have been spatially extended and 
of an extraordinary strength. 

If, on the other hand, the nucleus of the galaxy is in B2, 
as suggested by the near-IR colors, then this would imply that 
Bl harbors not only massive star forming regions but also a 
"runaway" super-massive black hole (SMBH). Did this SMBH 
originate from a companion which was accreted in the sys- 
tem at an earlier stage and it is now accreting gas becom- 
ing visible as an AGN? We do observe a high concentration 
of gas outside - but close to - the Bl and B2 regions of 
NGC 3690 and their kinematics imply that they are tidally con- 
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necte d with streamers originating from IC694 dCasoli et alJ 
1999) . lAlonso-Herrero et alJ (J2000) though derived a dynam- 
ical mass for Bl (M > 4 x 10 8 M^) similar to the one found in 
B2 (M ~ 6 x 10 8 M Q , from lShier etafl d 1 9961) 1 which would be 
hard to explain if both of them had always been part of a single 
galaxy. Our mid-IR observations do not allow us to discrimi- 
nate between the two interpretations. Only high spatial resolu- 
tion information on the kinematics of those components will 
help us shed some light to these issues. 

In any case, NGC 3690 has been the host of extraordi- 
nary star-forming events and exhibits a distinct mid-infrared 
spectral energy distribution unlike any other in other extra- 
galactic objec ts, even unambiguous A GNs such as N GC 1068 
jLeFloc'hetalJ200ll) or Centaurus A dMirabel etalJl 999). 

6. Conclusions 

Our mid-IR spectroscopy of Arp299 over the range 5 to 16 /urn 
has revealed new information on the properties of this intrigu- 
ing interacting system. Based on our analysis, we conclude 
that: 

1 . Large quantities of dust are present in the system as indi- 
cated by the hot continuum observed above 10 /im in all 
sources and by the high extinction derived by our model. 
We find that, assuming a uniform foreground screen model, 
we can reproduce the mid-IR SED and the 9.7 fim silicate 
absorption band with an A v of ~ 40 mag in the central 
region of IC694 (source A), ~ 10 mag for source C, and 
~ 50 mag for C', the dominant sources in the galactic disks 
overlap region. A best fit with an Ay of ~ 60 mag is de- 
rived in the central region of NGC 3690 (source Bl) using 
a homogeneous mixture of dust with radiating sources. 

2. The nuclear region of IC 694 harbors enshrouded massive 
star forming activity. No evidence for an active galactic nu- 
cleus is found in our mid-IR spectra. 

3. We confirm that C and C are sites of recent starburst events 
based on the elevated values of the [Nem] 15.6 /mi flux and 
the high [Nein]/[Nen] ratio, which is similar to those of 
classical young starbursts. 

4. Source B is characterized by a unique mid-IR spectrum 
which is dominated by a hot dust continuum extending 
shortward of 5 fim, and which is the strongest of the 
Arp 299 system. Given the ample presence of dust and the 
strength of the radiation field, this continuum is consistent 
with the recent Chandra/Beppo-SAX results which favor 
the presence of a Compton thick AGN at this location. 
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